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We compute the exact partition function for pure continuous Yang-Mills theory on 
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function is trivial, while in the strong coupling phase (large A) it is expressed in terms of 
elliptic integrals. We expand the strong coupling result in a double power series in 
and g^A and show that the terms are the weighted sums of branched coverings proposed 
by Gross and Taylor. The Wilson loop in the weak coupling phase does not show the 
simple area law. 

We discuss some consequences for higher dimensions. 
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It is an old and famous idea, that large QCD is equivalent to a string theory. 
Many approaches have been tried to make this precise. One of the most promising is to 
interpret the diagrams of the strong coupling expansion of the Euclidean lattice theory as 
string world sheets. This combines the great advantage of the finite strong coupling 
expansion, that confinement is already present at leading order, with the great advantage 
of large A^, that we have a free string theory at leading order. This expansion was proposed 
for any dimension in and elaborated in 

Early hopes for the use of the strong coupling expansion were dashed by the discovery 
of Gross and Witten and of Wadia [^0] , that even the simplest integrals involving the 
Wilson action, such as 

'dt/e^/«'*'-(^+^^)trt/ (1) 



(in two-dimensional Yang-Mills theory, this is a Wilson loop enclosing a single plaquette 
of area 1), were non-analytic in the coupling constant (see for the complete treatment 
of the problem). 

Such behavior is possible because we have 0{N) degrees of freedom, and these inte- 
grals are dominated by a saddle point. Just as for the infinite volume limit in statistical 
mechanics, several saddle points may exist, and the large limit picks the one of low- 
est action for a given value of the coupling. The consequence in higher dimensions was 
that integration formulas used in developing the strong coupling expansion were only valid 
down to a critical g^. The expansion had no validity in the physical regime of weak bare 
coupling. 

It was not then clear whether this transition is a fundamental barrier to combining 
strong coupling and large A^ or just a reason not to use the Wilson action in this context. 
Although we now know much more about large A^ phase transitions, after having studied 
them in depth for their application to two-dimensional gravity and string theory (this one 
was studied in 0), this point is still not clear. It is therefore interesting to modify the 
action and see if this modifies or eliminates the transition. The only obvious constraint is 
the very weak one that for small lattice spacing, the tr F'^ term should be present in the 
action, with higher dimension terms not unnaturally enhanced. 

The point was made by Migdal that in two dimensions, one could exactly integrate 



out a link common to two plaquettes, and compute an exact renormalization group trans- 
formation. The Wilson action is not a fixed point of this transformation; it and generic 

1 



nearby actions flow to the "heat kernel action," more properly described as a Boltzmann 
weight for a plaquette with holonomy U : 



ZHk{U;9'A) = J]dimi?e-^^^^^(^)/2^Xi?(t^) (2) 

R 

which is the heat kernel on the group manifold G{U, 1; At = g'^A/2N). This was also pro- 



posed in [|TT|] for its other nice properties, among them being that it gives exact equivalence 
between Euclidean and Hamiltonian lattice formulations, and that using the heat kernel 
action, the above Wilson loop expectation value suffers no large transition. From this 
one might conjecture that the transition was a lattice artifact, and that using the heat 
kernel action gives us continuum answers with no transition. (Of course this action, used 
for the plaquettes of a regular lattice, is not an RG fixed point in D > 2, so even if it 
worked for D = 2 the conjecture would not be proven.) 

Already in two dimensions the continuum QCD observables, Wilson loop averages: 

W{CuC2, Ck) =< n Pexp[i I dx^A^{x)] > (3) 

(where Cj are arbitrarily intersecting and selfintersecting contours on the infinite plane), 
show a great deal of structure. These were first computed by [jT^ for the U (N) gauge group 
with N CO, and generalized to any in [|T3[ and to the lattice version of the theory 



(see also [|T^,0). This calculation was based on the renormalized two-dimensional version 



of the Makeenko-Migdal loop equations established in and (in modern language) 
on the invariance of the results under area preserving diffeomorphisms. 

The most striking feature of the Wilson averages in the limit N —>■ oo, observed in 
was their "stringy" character: each of them was shown to be a sum over all possible 
surfaces of the minimal area (without folds) spanned on the contour ( "soap films" ) , of an 
exponent of minus area of the surface (area law) times some polynomial of the areas of 
domains forming this surface. The geometrical interpretation of the polynomials was not 
known at the time. It became clearer from the paper of I. Rostov (see Appendix A 
there) where it was demonstrated that they come from the statistics of surfaces (coverings 
with boundary) having branch points and cuts, connecting various sheets of a surface. For 
recent developments in this direction see . 



Another nontrivial quantity is the partition function on a two dimensional compact 
manifold of area A and genus G. A formula in terms of the sum over representations of 
the gauge group, 

Zoig'A) = 5^(dimi?)2-2Ge-.^AC.(i?)/2Ar_ 

R 



was found in [^. The sum for the case of the group U{N), say, goes over all Young 
tableaux characterized by the components of the highest weight {tii, 712, un} which are 
the integers obeying the inequality: 

00 > ni > ... > njsf > —00 (5) 

So we see that it is still a complicated multiple sum which takes some effort to calculate 
in particular cases. 

A nice interpretation of this partition functions in terms of minimal coverings, similar 



to those of [|T2[ and , was given by Gross and Taylor. |21|,^ It was noticed that the 
partition function can be written in terms of a sum over minimal coverings of a manifold 
of a genus G and area A, where different sheets of a covering are glued together with 
elements called branch points, tubes and omega-points. Using an inequality known by 
mathematicians, about the possibility of minimal covering of the manifold of a genus G by 
a surface of a genus g. Gross found that many terms of the topological expansion in 
the free energy F{A, N) = log Z{A, N) are equal to zero. So, for G = 1 he found that 
the 0{N'^) (spherical world-sheet) contribution to F{A,N) is zero. However, for the next 
order (torus topology of coverings) the sum over minimal coverings is infinite and can be 
given in terms of the Dedekind function: 

F = -2 logr^iiA/A-K) = - 2 J] log(l - e^^/^) (6) 

n>l 

(the constant is of course a choice of ground state energy.) 

For the spherical topology of the space G = 0, the result is more complicated to 
obtain. Its geometrical interpretation contains all of the additional elements mentioned 
above, and a purely geometrical derivation looks tricky. 

In this paper we present the explicit result for the leading order (planar) contribution 
to the free energy of the Yang-Mills theory on the two-dimensional sphere. For large 
area of the sphere (compared to it nicely fits the interpretation in terms of minimal 



coverings, down to the phase transition point g'^A^ru = tt^, where the sum over coverings 
is divergent. In the phase of smaU g'^A the result is trivial. 
The partition function is 

Zg=o(^, N = oo)= exp[iV2F(A)] = ^{dimRfe'^^^^^'^ (7) 

R 

where we measure the area A in units of l/s*^, and for the group U{N) we have 

AT 

C2(i?) = ^n,(n, -2z + Ar + l) 

(8) 

dimi? = rr(i-^^^^) 

i — j 

and the sum over the representations R has to be understood as a multiple sum over 
integer variables ni, tiat obeying the inequality (^. f Now, in the large limit, nothing 
prevents us from using continuum variables: 

obeying now the inequality: 

n{x) > n{y), if x <y. (10) 

It is convenient to change the variable to 

h(x) = -n{x) + X - 1/2 (11) 

and to write formally the following functional integral representation for the partition 
function: 

ZoiA)= J Dhix)exp-N^Seff[h{x)] (12) 



where 

SeffiHx)] 



[ dx [ dylog\h{x)-h{y)\ + ^ [ dxh'^{x) ~ A/24: (13) 
Jo Jo 2 Jq 



f Strictly speaking these are only representations whose U{1) charge is correlated with the 
charge under the center of SU{N), in other words of SU{N) x U{1)/'^n ■ Since a single represen- 
tation will dominate in the final answer, it is also correct for U{N). 
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and h{x) obeys, according to and (^TJ) , the inequality 

> 1. (14) 

x-y 

One of the main observations of this paper is the need to respect this condition in calculat- 



ing (|T^) , which will lead to nontrivial consequences such as a new large N phase transition 
and the existence of the large area (strong coupling) phase. If we introduce the density of 
the boxes in the Young tableau in terms of the variable h 

uih) ^ ^ (15) 

with the normalization 

dh u{h) = 1 (16) 



the condition (|14D can be simply rewritten as 

u{h) < 1, for any h (17) 



Since we have a large parameter in front of the S^f j in (|T2D we can try to apply 
the saddle point approximation, which means that we have to solve the equation on h{x) 

6S,ff[h{x)] ^ ^ ^^g^ 
dh{x) 

Let us ignore for a moment the constraint (^7]) . This will lead us immediately to the 
integral equation 

^h = p[p^ (19) 
2 J h-s ^ ' 

which is precisely the same as for the distribution of the eigenvalues in the hermitean 

gaussian matrix model. The solution is the well-known semi-circle law of Wigner: 



From here we obtain immediately for the derivative of the free energy with respect to 
the area A 

dS^ffm tr 1 _ 1 1 

^ dA- -< 2iV^ > "24 - 2l " 24 ^^^^ 

or 

F{A) = ^^A-\\ogA (22) 



This result is already clear from (^2]) where the log A term can be obtained by simple re- 
scaling of the continuous field h by \/A, after which the gaussian integral will not depend on 
A. This result was obtained by the same reasoning in |^^, where the inequality (^7|) was 
ignored, and hence, the existence of the second, strong coupling phase was not noticed. 
Let us explain this result from another point of view. Clearly from (|^), we have 

ZG=oiA) = ZhkiU =1; A)= G{U = 1, C7' = 1; A/2N). (23) 

For small A, we really are interested in the small time asymptotics of the heat kernel on 
the group manifold. These are well known for arbitrary group manifold H: 

^""-'^^^ ^ [£a) +0{e-'-^'^) (24) 

where the correction is the sum of over all closed orbits on H. 

We see the origin of the term — | log A, and that the corrections are suppressed as 
^-N/A^ hence vanishing in the large limit. 

Remembering now that we have the constraint ([T^), we note that the trivial solution 
( pOD is possible only for small areas A: 

A < A,r^t = TT^ (25) 



What happens for A > Acrit! We still have to solve the saddle point equation 
but in the presence of the boundary conditions (|5|) , or, which is the same, (|T^) . The 
only thing which may happen is that a finite fraction of the highest weight components 
ni, ...TT-AT will condense at the boundary of the inequality, namely 

^fc+i = = •■■ = J^AT-fc = (26) 

whereas all others are non-zero (see fig. lb and compare it to fig. la, where a typical Young 
tableau for the weak coupling phase is presented). 

On the language of the density for the continuous variable h it means that 

uih) = 1, for -b<h<b 

(27) 

= u{x) elsewhere 

where 6=1/2 — k/N is finite in the large N limit and u is some nontrivial function to be 
found later. 
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In fig. 2 the Young tableau is presented in both phases for the variable h{x). We are 
trying to find here a so called two-cut solution, similar to those observed in the hermitean 
matrix models with double well potentials [p4 |. 



Let us substitute the ansatz ( pTf ) into (|T^) . We get the following equation on u{h): 

As we see now, the condensate of the zero highest weight components induces an extra 
logarithmic term in the equation. In the language of the equivalent hermitean one matrix 
model it would mean that we have the effective matrix potential whose derivative is the 
l.h.s. of ( PBD . This potential clearly has two wells separated from each other by the cut 
of logarithm. The eigenvalues fill these wells always to the top. All the eigenvalues which 
spill over the top (say, by increasing A) form the condensate. 

We introduce, as usual, the function of the complex variable h 

f{h) = I ds^^ (29) 

whose imaginary part is nu{h). The solution of the corresponding two-cut Cauchy problem 
p5[1 is given as a contour integral 

1 , r ^ As - log ^ 

f{h) = — -v/(a2-/i2)(62_/,2) i^^ ^ ^h+b (^3Q^ 

where the contour of integration encircles the cuts of the square root but leaves aside the 
singularities of the nominator and of the pole at s = h. The limits a and b have to be 
found later from the condition of the correct behaviour of f{h) for h oo. 

Now, by inflating the contour we catch, instead of the cuts of the square root, the 
singularity at the pole and the cut of the logarithm. This gives 



Let us note that the imaginary part of the logarithm in the r.h.s. of (|31] ) is exactly 
equal to the minus tt times density of condesate of h's in the interval [—6, b]. So it is clear 



that the last term in (^) represents the full function u{h) defined by (^Tj) . The latter is 



expressible in terms of the complete elliptic integral of the third kind 11 [^, x] 



„rM - 1^ l {a + h){b + h) 2b h-a 2v^ 

^ ^" 7r6 + aV {a-h){b-hy^^a + bh + b' a + b^ ^'^ ' 



In order to find the parameters a,b and the whole free energy, it is better to use the 
asymptotics for the large h of ( pTD 

A 1 

f(h)=h(-- ds , ) 

g2 _ g^+b^ 



^'^-'(''^IJ'^/W^^W^)^ (33) 

+ h-^2b+ I ds , 2 ^^)+o(/,-5) 



and compare it with that which follows from the definition (|2^) : 



f{h) = • /i - (1 - 2b)h-^ - {F'{A) + l/24)/i-3 + 0{h-^) (34) 

This comparison gives (the elliptic integrals are in the appendix) the following results: 
The first derivative of the free energy in the strong coupling phase, expressed in terms of 
elliptic integrals with modulus k = b/a, is 

F'{A) = -a" - —a^k'^ - — + —a^k'^A (35) 
^ ^ 6 12 24 96 ^ ' 

where the modulus is related to the area by 

^A = {2E - k'^K)K, (36) 
the complementary modulus k'"^ = 1 — k"^, and 

a = 4K/A. (37) 

This solution represents the strong coupling phase of our theory, namely for the area 
of sphere A > n"^. It is easy to check that at the point of transition Acrit = t^"^ the two 
solutions coincide completely, even for the distribution u{h) of boxes in the Young tableau. 
Let us calculate the order of this transition. 

Series expansions become easier in terms of theta constants for a torus of complex 
modulus r. The equation (|3^ ) becomes (the relevant identities are in the appendix) 

A = 8EK - Ak'^K^ 

,2 



^{el{0\r) + el{0\r) + 2E,{r)) (38) 
n^{l + 8q - 8q^ + 32q^ + . . .) 
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where q = e*'^'^, and the critical point is the hmit r — zoo. 
Now 

^0 as r ^ zoo 

so we see that the transition is higher order. Inverting ( |38|) and substituting, 



1 ^A-A,^' 



+ 



Thus the phase transition is of the third order, hke the weU known Gross- Witten- 
Wadia phase transition for the lattice two dimensional multicolour gauge theory. However, 
in spite of some similarities of these two transitions, the one found in this paper happens 
already in the continuum version of the theory, so we cannot say that it is a lattice artifact. 

The transition also bears some similarity with the Berezinski-Kosterlitz-Thouless tran- 
sition of condensation of vortices on the world sheet of one-dimensional string theory com- 



pactified on a circle . In the language of the corresponding matrix quantum mechanics 
the point of the phase transition also corresponds there to the disappearance of the gap in 
the characteristic Young tableau for the U{N) representations of angular matrix variables 



27 



We can make contact with the results of Gross and Taylor by expanding the answer 
about g'^A = oo. Although this is a singular point, the form of the singularity allows a 
well-defined double expansion in e~"^/^ and Ae~"^/^, as will emerge in the following. From 
21j], a reason to think that the expansion is unambiguous is that each term e""^"^/-^ has 



a coefficient polynomial in A and of order 2n, a property we would certainly lose if we 
expanded the exponentials in some other way. A better argument requires knowing the 
analytic structure of F{A), to which we turn. The strong coupling limit was r — 0; since 
series expansions of the theta constants are in e*'^'^, clearly we want to make a modular 
transformation. This can be done by taking K ^ K' and k ^ k' and then going to theta 
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functions with modulus r — ^ zoo in the strong couphng hmit. Using Legendre's relation 
for E' gives 

= -k^K'"^ + 2K'/K{tx/2 + KK' - EK') 

= -^ir + ^^(^3'(0|r) + elm - '2E2{t)) ^^^^ 
d 

= -TxiT - 2niT^— log^4(0|2r) 
or 



-TxiT + (27rzr)^i?(27rzr) 



Now 



R{2mT) = 2 5^ - 



g4ni7rT 



n>l 



(42) 



so in the limit A ^ oo we have r = iA/ A-k with corrections exponentially small in A. To 
get a double expansion in exp— A/2 and Aexp— A/2, we express everything in terms of 
theta functions with modulus 2r, and solve for 27rzr in terms of A and the exponentially 
small R: 

2mT = ^(1 - Vl + 4RA) 

^ -^A s{AR) (43) 

= -^A{1 -RA + 2R^A^ - hR^A^ + lAR'^A'^ + . . .) 

We would then successively substitute A for r in R. 
Rewriting F'[A) in the same way, we find 

F'{A) = + Y^{AR)\et{Q\r) - \9l{Q\T)) + ^s{AR)'91{Q\t)A 
= + ^^s{AR)\9t{Q\r) + el{Q\r)) + ^s{ARYel{Q\r)A 

= + s{ARfFl^{2TiiT) + A s(Ai?)^F^(27rir). 

We see that the analytic structure of F' in terms of w = e~^l'^ is not so simple; 
however the branch cut in (^) is away from the origin, and near the origin we have a sum 
of terms {Vo'gw)'^ j raiw) with each analytic. If we were only given the function F{w), 
we could isolate these terms by combining its values on the sheets F{e'^'^^^w); at each order 
in w only finitely many fm contribute. This would fix the double expansion uniquely. 
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Using d?^) and (^) , we can get the terms at a given order in the double expansion 
to aU orders in e""^/^ by taking AR smaU; then (here prime is always d/dA = —4q d/dq) 

F'{A) = - ^ + Fl,{-A/2) + A{F[ - 2i?F^) 
+ A^ {-4RF{ + hR^Fl^ - RFl^) 

~ 2-^ I _ g-(2n-l)^/2 

(45) 

(1, >^ ne-"-^/^(l - (-e--^/^)") . >^ (2m - i)e-(2"^-i)-4/2 

"'" Ig^ '^ Z^ 1 - e-"^ ^ 1 _ g-(2m-l)A 

^ n>l m>l 

me-"^-^/^ , 1 ^ (2n-l)e-('"-^)-^/^ ' 

1 _ g-m^ ^ ]^ _ g-(2n-l)A/2 

m>l n>l ' 

+ 0(A2). 

We have compared this with a direct expansion of the formula (|^) to 0(exp(— 2A)): 

F{A) = 2e-^/2 + (-1 - 2A + -A2)e-^ + (^ + 4^2 _ 8^3 ^ \j^^y-zA/2 ^ _ _ _ 

2 3 3 3 

-F'{A) = + (1 - 3A + -A^)e-^ + (4 - + 2^^ - -A^ + 1^^)6-3^/2 + . . . , 

2 3 2 

(46) 

and with the 0{A^) and 0{A) terms to much higher order, and found complete agreement. 



One can also reproduce this result by the direct expansion of (|35|) . f 

We also checked this expansion from Gross and Taylor's rules, dropping the terms 
involving "tubes" and "handles" and proportional to powers of A (as appropriate for 
U{N)). 

One technical conclusion we can draw from the solution is that and A are not 
really the natural expansion parameters in the problem. One way to think of this is in 
terms of the formalism of [1^]. The perturbation is a higher derivative operator, and it is 
perhaps surprising that this even has a non-zero radius of convergence. Evidently it does, 
and much of the effect of the perturbation can be expressed as a "renormalization" of the 
modulus of the cylinder from A to the variable r determined as above. 



f We thank J.-M.Daul, who reproduced this expansion starting from the integral representation 
|3^ ) up to 0(exp(— |A)) and found complete agreement. 
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In Gross and Taylor's language the power-like expansion in A is given by inserting 
branch points on the string world- sheet. Evidently this expansion diverges at Ac- 

We might hope to see the transition from the weak coupling side by noticing that 
some expectation value obtains an impossible value in the continuation beyond A^ just 
as for the GWW transition, the continuation of the strong coupling spectral density to 
small g was no longer positive. The easy observables to compute here are expectation 
values of the local operators tr i?^ (E is the electric field), given by J dhu{h)h^ . Another 
possible way to see the transition would be to see the sum of terms e~^^^ in (|2j) diverge 
for sufficiently large A. 

Another quantity revealing of the weak coupling phase is the Wilson loop, separating 
regions of area Ai and A2. This is again simple in terms of the heat kernel: 

Wg=o{Ai,A2)= j dUG{l, U; Ai/2N) ^tr {U + U+) G{U, 1; A2/2N). (47) 

Prom (|2^) we learned that in the weak coupling phase, only the leading classical solution 
contributes in the heat kernel. Thus we can take the expression of |]TI| and drop the 
winding terms: 

G(l, U; A/2N) = ATJ] -4^7^\t^"^''^^^ (48) 
where e*^* are the eigenvalues of U. Since the invariant measure is 

[dU = jUde, JJsin^ 1(0, - e,), (49) 



( ^71) becomes 

Wg=o{A^,A2) = jUde, 

i<j i 

which is again an expectation value at a semicircular saddle point: 

Wg=o{A^,A2)= (^-^^ Ji(v^) 
= 1-X/8 + ... 



(51) 



where x = ^1^2/(^1 + A2). 

This is perhaps a peculiar result. It is positive for all Ai + A2 < Acrit but becomes 
oscillatory for (unphysical) large A. Even stranger, its asymptotic behavior is cos{^/A)/A. 
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It would be very interesting to calculate the Wilson loop in the strong coupling phase. 
Knowing the result for the simple loop, the two-dimensional renormalized loop equations of 
||12|| would then determine all loop averages. It is clear that these will obey the same Gross 
and Taylor rules of the large Ai, A2 expansion (but the coverings will have the topology 
of the disc now) . In the limit of large total area of the sphere we would reproduce all the 
results of for the Wilson average on the infinite plane. Only the coverings which do 
not wind over the sphere will survive, and they are precisely those which were observed in 

Returning to the general consequences of the phase transition, the conclusion for the 
relation of the strong coupling expansion to string theory is that string rules derived from 
the heat kernel action by expanding about = 00 (in ||2^, the expansion was in exp —g'^), 
do not give the correct answer on a small two-sphere. Now if we were two-dimensional, 
we might not care about this - the string rules DO give the right answer in the A — ^ 00 
limit, and all we seem to need for this is that the overall area of our universe be large. 
Expectation values for a Wilson loop enclosing an arbitrary area in this universe will be 
given correctly. Certainly the 't Hooft model of mesons in QCD2 has no large transition 
in infinite Minkowski space. 

For four dimensional physicists, however, this transition looks like a real problem. The 
only precise way we know to define the strong coupling expansion is to start on a lattice, 
and take the continuum limit as defined by Wilson. For QCD this will require taking the 
bare coupling to zero in the way prescribed by the RG, so our expansion must make sense 
at weak coupling. We therefore need to choose an action for which there is (among other 
constraints) no large transition. Although we do not know if the transition we find 
for the heat kernel action persists in D > 2, the simple fact that any higher dimensional 
lattice contains embedded topological two-spheres puts the burden of proof on the other 
side - to show that somehow cancellations between terms in the higher dimensional series 
eliminate the transition. 



We thank T.Banks, E.Brezin, D.Boulatov, J.-M.Daul, J.Distler, P.Ginsparg, D. Gross, 
C.Itzykson, I.Kostov, H.Neuberger, S.Shenker and M.Staudacher for enjoyable discussions. 
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Appendix A. Elliptic Integrals and Theta Functions 

The basic integrals required are 

dX 2 



K{h/a) (A.l) 



v/(a2-A2)(62_A2) a 



and 



-==_=== = -a[{2a' + b')K{b/a) - 2{a' + b')Eib/a)] (A.3) 

in terms of the standard complete elliptic integrals with modulus k = b/a (e.g. as in ||29|| ) 
We will write the results in terms of the complementary modulus k' satisfying k'"^ = 
1 — /c2, using K{k) = K'{k'), E{k) = E'{k') and re-express them as theta constants, for 
the strong coupling expansion. From (13.20) we have 



A: = ^2(0|r)V^3(0|r)2 

k' = ^4(0|r)V^3(0|r)2 

K{k) = lel{0\r) 

2 (A.4) 

K'{t] = ^«|(0|T) 



36'4(0|r) ' ' 12K{k) e[{0\T) 



where q = e*"^"^ = exp{—TTK'/K). 

Some standard identities involving these functions: 

KE' + K'E - KK' = -TV 

2 

^2 d 

-Or, = 4Tci-—9„ Vn 



9z^2 Qt- 

e[{0\T) = 7r^2(O|r)03(O|r)^4(O|r) = 27ir]{Tf (A.5) 
^3^(0|r)=^2^0|r) + ^|(0|r) 

-jrf^ = 4m-log9[{0\r) = 127rz-log,7M = -n'E^ir) 
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where E2 is the normahzed Eisenstein series. We will need as well 



-Kidr "^4(0|r) ^^g^ 



TTi dr ^4(0|r) 

proven by checking that both sides are modular forms of weight 2 and level 2, and low 
orders in the g-expansion. 

All of the relevant theta constants can be expressed in terms of ones with modulus 2r: 

elm = {91{0\2t) + 9l{0\2T)r 

et{0\T) = {el{0\2r) - el{0\2r)f 
.4,^,.^ , .4,.,_^ 8 d ,_^2(0|2r)^3(0|2r) 



^3*(0|r) + ^|(0|r) = - — log 



nidr ^ 0j{0\2T) (A.7) 

12 d , ^4(0|2r) 
TTT log 



ni dr ??(2r) 

E2{t) = -—log04{O\2T)Tj{2T) 

m or 

The first two follow from 13.23.15; the third uses these, ( [A.6| ) and Jacobi's identity 
(line 3 in ( [A.5| )); the fourth follows from substituting the product representations below. 

We will then need series expansions of these and their r-derivatives. These are best 
derived from the logarithmic derivatives. (Here q = e^'^'^): 



1 d , ^4(0|2r) 1 {2n~l)q 



4Tvi dr ^ ?7(2r) 24 ^ 1 - Q4n-2 



4n-2 



n>l 



n>l ^ 



,4n-2 



n>l ^ 



m>l 

.Am 



E2{2t) = 1 - 24 5^ - 



m>l ^ 



Finally, we use product representations like 

^4(0|t)= n(i-9''^)(i-9'""')' 



n>l 

n>l 



(A.9) 
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